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Abstract 



The results of numerical calculations of linear polarization from magnetized spherical optically 
thick and optically thin envelopes are presented. We give the methods how to distinguish magnetized 
optically thin envelopes from optically thick ones using observed spectral distributions of the polar- 
^v^ , ization degree and the positional angle. The results of numerical calculations are used for analysis of 

\_C) • polarimetric observations of OB and WR stars, X-ray binaries with black hole candidates (Cyg X-1, 

f^ ' SS 433) and supernovae. The developed method allows to estimate magnetic field strength for the 

^^ , objects mentioned above. 

o 

-g^: 1 INTRODUCTION 

O I The existence of magnetic fields in hot stellar atmospheres and envelopes as well as in accretion envelopes 

-4— > ' of quasars and AGNs gives rise to the Faraday rotation of the polarization plane of the radiation from 

jrt ! these objects. The angle of the Faraday rotation ip may be vifrittcn in the form [1]: 



iP {n, B) = -dTT cos e, (1) 



C^ ■ S=^■-■^=0.8X^'i^lm)B{G). (2) 
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where n is the direction of wave propagation, O is the angle between n and magnetic field B, tt is 
the Thomson optical length [tt = N^cttI), N^ is the electron number density, / is the geometric path, 
aT — 87rrg/3 is Thomson's cross-section, r^ — e^/rac^ k, 2.82 x lO^^^cm is the classic radius of electron, 
ljJb — \e\B/7nc is the cyclotron frequency, uj = 2t:c/X = kc is the frequency of radiation (ljb/lo « 
0.93 X lQ'^\{iJLm)B{G)). 

At large distances from a star or a quasar magnetic field is approximately dipolar because other 
multipoles fall off with distance more rapidly than the dipole component. Because of small value of 
Thomson's cross-section ctt the optically thick envelopes seem be seldom the case. The expanding matter 
just after the supernova explosion, when the concentration of particles is yet large enough, is possibly a 
case of such optically thick envelope. Another example may be the accretion structures near quasars and 
AGNs. 

Magnetospheres around the black holes of different masses have to play a very important role in 
formation of various structures of outflowing matter (winds, expanding coronae, jets, etc.). Possible 
existence of such magnetospheres was confirmed recently both by spectroscopic (Robertson and Leiter 
[2]) and by polarimetric (Gnedin at al. [3]) observations of X-ray binary systems with black holes. 

Now the Blandford and Znajek mechanism is commonly accepted which provides the energy extraction 
from a rotating black hole by means of magnetic field [4] . 



Later this idea was developed in papers [5-7]. In these papers the model was considered where the 
toroidal electric current at the inner boundary of the accretion disk generates the dipole magnetic field 
connecting together the black hole horizon and the region of the accreting disk beyond its inner boundary. 

The plasma outflow from the magnetosphere and the inner part of accretion disk gives rise to an 
envelope. The radiation scattered in the envelope acquires the linear polarization. The region of the 
formation of broad emission lines observed in AGNs is a case of such envelope. 

In the absence of magnetic field the scattered radiation has to be non-polarized provided the scattering 
matter in the spherically shape envelope is symmetric relative to the line of sight. If the magnetic field 
of the object is not symmetric relative to this line, the outgoing radiation (singly or multiply scattered) 
has to undergo influence of the Faraday's rotation effect and becomes to be polarized. The integral 
polarization of radiation is equal to zero if the dipole magnetic momentum M is coaxial with the line of 
sight or in the case of B = 0. 

2 OPTICALLY THIN ENVELOPE 

The integral polarization scattered in an envelope is simplest to calculate if the envelope is optically thin. 
In this case we can take into account only single scattered radiation and the Faraday rotation during the 
followed propagation to an observer. This was made for various types of envelopes in papers [8-10] , where, 
however, was accepted that a star is a point-like source of non-polarized radiation. This assumption gives 
rise to some overestimation of the polarization degree (more detailed, see below). Maximal disturbance 
of axial symmetry of the Faraday's rotation picture occurs when magnetic dipole M is perpendicular to 
line of sight n {{)m = 90°). In this case the integral polarization from the whole envelope acquires its 
maximum value and the polarization plane coincides with the plane (nM). 

The polarization Stokes parameters Fq and Fjj, according to [1, 8, 9], are described by the following 
expressions: 



L 3 _ f,,,N,{r),L m^... 



^^("^^"i^'Te^^^/ '^^^V'"7^'^"^'°'^^^ + ^^' ^^^ 



F^(„) = --^ • yI^-t/ rfF^^l- ^sin^^ sin2(^ + ^), (4) 

where L{erg s~^ Hz^^) is the luminosity of the central star, R is the distance from the star to the 
telescope, the angles i? and ^ determine the direction of the radius vector r(r, t?, cp); Rg is the radius of 
the star. Integration yields over the visible part of the envelope volume. In (3) and (4) we accepted that 
X-axis lies in the plane (nM). 

Near the stellar surface the polarization of the scattered radiation is small because of the non-polarized 
star's radiation is falling to the volume element dV practically isotropic from the solid angle « 27r. In the 
book [1] (chapter 4) the formulae are given which describes this effect quantitatively. There was considered 
the commonly used model when the intensity /(/i) of outgoing radiation is approximated by the formula 
/(/i) = A + Cji. Remember that /i is the cosine of the angle between the normal to the atmosphere's 
surface and the direction of the escaping photons. In this paper we consider only the isotropic case 
(C = 0, Lambert's law of radiation) when the effect of depolarization of scattered radiation near the 
stars's surface is maximal (most profound). This gives rise to the factor -^/l — _R|/r^ in formulae (3) and 
(4) as compared with the analogous formulae in [8-10] where a star was considered as point-like source 
of non-polarized radiation. The non point-like character of the radiation source was taken into account 
first in [11]. 

We restrict ourselves by two types of the electron number density A^e(r) in an envelope: N^ir) = 
Nq — const and A^e(r) = iVo(-Rs/r^). The explicit formulae for Faraday's rotation angle ip have the form 
[8]: 

For the envelope with Ne{r) ~ const 



O T O / O O 

^ = r>/ ^ ^T\ o jcost^m [cos^ -Jl -sin^^] + sin??„ cos(^sin^(l -)}. (5) 

Here, rj = Rq/Rs is the ratio of the envelope radius Rq to the radius of a star R^; p — r/Rg, Sg = 
0.8A^M/i?g is the Faraday's depolarization parameter (2) at the magnetic equator of the star, M is 



Table 1: The coefficients Kq, Ku and the maximal polarization degree for the various values of the 
parameter i] (the model of a non-point like source of radiation). 



IOOKq 0.747 1.116 1.362 1.376 0.830 0.325 0.078 0.028 0.013 

IOOKtt -2.644 -2.944 -2.593 -1.638 -0.375 -0.048 -0.002 -0.000 -3 10"^ 






Table 2: The same values for the model of a point like star. 



V 


1.02 


1.05 


1.1 


1.2 


1.5 


2 


3 


4 


5 


IOOKq 
100 Ku 


8.018 
-28.32 

8.136 
1.30 


7.593 
-19.67 

8.728 
1.35 


6.631 
-12.27 

9.310 
1.40 


4.892 
-5.592 

9.830 
1.55 


2.065 
-0.869 

10.03 
1.90 


0.662 
-0.086 

9.241 
2.55 


0.139 
-0.003 

7.403 
3.65 


0.048 
-0.000 

6.609 
6.30 


0-022 
-0.000 

6.234 




^^. -e.. 


10-00 



magnetic dipole moment of the star, Tenv = N()aT{Ro — Rs) is Thomson's optical depth of the envelope, 
dm is the angle between M and line of sight n. 
For the envelope with Ne — No{Rl/r^): 

■4) = '^^'"1 {cos^m [-(1 - cos??) + ^(1 - cos^??) - -(1 - cos'"^^)] + ^ sin t?„ COS (^ sin^ t?}. (6) 
2/?'' sin 1? 3 5 5 

Here, Tenv = NqGt Rs IS the optical depth of the envelope. Remind that for single scattered radiation 
{Tenv '^ 1) we neglect the terms ~ r^ and higher. On this approximation when calculating the polarization 
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degree pi = \/Fq + F^/Fj one can take for the flux Fj the value L/AttR^ 
It follows from (3) and (4) that Fq and Fjj have the form: 

FqM = FiTenv fQ^ui^sTenv , l^m, V) ■ (7) 

For the envelope with N^ ^ r^^ the variable r/ in (7) is absent. 

For weak magnetic fields or large wavelengthes, when -0 <C 1, the fq^jj functions have the following 
asymptotic form: 

/q C:^ Kq{Ss Tenvf sin^ ??„, 
fu ^ Ku ((5s TenvY COSdm STC? l?„. (8) 

The case ■0^1 corresponds to SsTenv ^ 1. It means that \fu\ "C |/q|, i.e. the wave electric field 
oscillates near the plane (nM) and the polarization degree pi « Tenv\fQ\- The asymptotic formulae (8) 
take place up to SsTenv ~ 1- 

The numerical coeflicients Kq and Ku for the case Ne ~ r^"^ are equal to 0.0006641 and -0.0000796, 
correspondingly. For the model of point-like source of radiation these values are considerably larger: 
0.0016136 and -0.0001873. In the case N^ = const the coeflicients Kq and Ku depend on the parameter 
rj = Rq/Rs- The numerical values of KQ{ri) and Kuifj), as functions of ?], for non point-like star, are 
presented in Table 1. It should be noted that for t] < 1.15(ii'Q(1.15) = 0.014136, i4:c/(l. 15) = -0.02088) 
the polarization of radiation decreases. This fact reflects the effect of depolarization of scattered radiation 
near the star surface, mentioned above. On the contrary the values of _fi'Q(?7) and Kuii]) for the point- like 
star model (see, Table 2) increase monotonously with 77 ^ 1. The negative sign of Ku denotes that the 
positional angle x (tan 2% — ful fo) is also negative, i.e. corresponds to left-hand rotation of the wave 
electric flelds from the plane (nM), for the line of sight n. 



The values pi/renv — \KIq + fu) (i"^ '^) ^"^^ ^^^ positional angle x (i'^ degrees), as the functions of 



^6sTenv, are presented in Fig.l. Remember that the X-axis lies in the plane (nM). The upper panels 
correspond to N^{r) ~ No{Rs/r)'^ and the lower ones correspond to envelope with Ne — const and 
77 = 5. All curves represent the model of non-point- like star. For the polarization degree to be obtained 
the right-hand curves must be multiplies by the optical depth Tenv < 1- 

For optically thin envelopes the maximum polarization pimax (at -dm — 90°) can acquire rather 
considerable values ^ 6Te„„(%) for the non-point-likc star (see Table 1) and ~ lOTenv {%) for point- like 
model of the radiation source (see Table 2). In these tables we present the values of V^sTenv corresponding 
to pimax values. 

It is notable that pi max is almost the same for various geometric depths of the envelopes (77 ~ 1.5 -^ 5 
for non-point-like stars) and (77 ~ 1.02-^3 for point-like model). Only the corresponding values of y/SsT^nv 
parameter change its position. Considering the fixed value r^nv for all envelopes we can interpret the 
mentioned fact as for every envelope there exists the effective layer with ip m 1 which gives the main 
contribution to the integral polarization. The interval of \/6sTenv, corresponding to pi max , increases with 
the increasing of the envelope radius _Ro- For 77 = 2 this interval (the width of the polarization spectrum) 
corresponds to ^/SsTenv ~ 2 -^ 5. For rj — 5 this width corresponds to y/Ssr^nv ~ 4 -^ 20. 

With increasing magnetic field this effective layer shifts from the stellar surface up to the outer 
boundary of the envelope. For further increasing of B the polarization is determined by the thin boundary 
layer where -0 < 1. This condition takes place at the optical depth of this layer r ^ I/Sq oc A~^, where 
60 = 0.8X^M/Rq is the parameter (2) at the outer boundary of the envelope. This estimation assumes 
that the total inner volume docs not give a contribution to integral polarization, i.e. SsTenv/v^ 3> 1. 
The intensity of the radiation scattered in this thin layer is proportional to t. It means that the integral 
polarization is also inversely proportional to Sq, i.e. pi ~ C(i?„i)?7'^/A^. The analytical formula for 
dm ^ 90° gives: 



l6Ss V V^ -^^ ' V^ 

Asymptotic dependence pi ex A^^ takes place for 5sTenv/'rf ~ 50 -f 100. The expression (9) was 
obtained without the exclusion of invisible part of the envelope, i.e. it gives rather overestimated value 
for pi. So, for ?7 = 2 the overestimation is equal to 11% at \/5sTenv — 1 and equal to 3% for ^/SsTenv = 10. 
For more vast envelope with 7] ~ b the corresponding values are 6% and 0.6%. For the model with 
Ne OC r^^ these values of overestimation are equal to 15% and 2%, correspondingly. If 77 ^ 1, the 
invisible part of the envelope are practically the half of all the envelope volume. Therefore, for this case 
we need to diminish the value (9) to 2 times. 

As it was showed in[l], a star can be considered as a point-like source of non-polarized radiation 
from the distances (1 -=- 2)Rs from the star surface. Therefore, the asymptotic values of polarization 
at SsTenv S> 1 are the same for both models of point-like and non-point-like sources, if, of course, the 
geometric depth of an envelope is larger than (2 -h S)Rs- 

For envelopes with N^ ^ r^^ the peak of polarization (at -dm — 90°) occurs at ^/dsTenv ~ 4.65 
and equal to 2.886Te„„(%) for the case of non point-like model. The point-like model for the same 
value \/6sTenv ~ 4.65 gives 3.878Tent)(%)- The peak polarization for point- like star model occurs at 
VSs'Tenv ~ 3.05 and is equal to 4.722x6™, (%). 

Asymptotic behavior of polarization at ^/SsTenv S> 1 and Ne ~ r"^ is determined also by outer layer 
of the envelope (r > r*) where V'(''*) = 0.55(TH.)T(rH.) < 1. The optical depth of this layer is equal to 
T{r^) — tenvRs/'''*! where 5(r*) is the value of parameter (2) at this boundary r = r,. Taking into account 
that Sir^,) ex 1/r^, we obtain r, ex vA. This means pi oc 1/vA. To calculate the polarization at large 
values of the parameter ^/SsTenv we can, as in previous case, take the numerical value at \/5sTenv = 10 
from the Fig.l and extrapolate it, using asymptotic dependence on A. For 'dm — 90° and Ne ~ r~^ one 
yields the formula: 

oc — =, OsTenv > 1- (10) 



'^ {SsTenvY/*' VA ' 

Thus, for Ne ^ r^^ the polarization degree decreases as ~ I/a/A- Using qualitative derivation of 
Pi{X) at \/SsTenv S> 1, givcu abovc, one can derive for Ne oc r~'^ the following asymptotic A-dependence 

PiOcA-2(--l)/(-+2). 



The polarization at dsTenv S> 1 is determined by the contribution of the envelope's volumes far from 
the central source. It means that both the polarization degree and the positional angle give an information 
about number density Ni,{r) and magnetic field B(r) just in this part of the envelope. So, polarimetric 
observations can determine (or to prove some accepted model) the distribution of free electrons in the 
envelope. 

3 OPTICALLY THICK ENVELOPE 

To calculate the integral polarization of radiation from optically thick envelopes one has to know the 
intensity / and the Stokes parameters Q and U of the radiation outgoing from an element of the surface 
of the envelope. We assume that the sources of radiation with the considered wavelength A are distributed 
far from the outer surface of the envelope, i.e. the parameters /, Q, U are the solution of the classic Milne's 
problem for semiinfinite plane-parallel atmosphere with taking into account the Faraday rotation of the 
polarization plane. This problem was solved numerically (see, [12, 13]) only for the magnetic field directed 
along the normal to the atmosphere. For the arbitrary direction of magnetic field the problem is too 
complicated and is not solved yet. Considering the dipole, i.e. non- homogeneous, magnetic field we are 
to know the solution of the Milne problem for arbitrary magnetic field direction. 

In the paper [14] the simple asymptotic formulae (at Sq > 1) for /, Q, U are derived for a number of 
radiative transfer problems, including the Milne problem. It is important that these analytic formulae 
are suitable for the arbitrary direction of magnetic field. They give exact solutions for large values of 
the depolarization parameter S (see (2)). These formulae correspond to the approximation when the 
intensity of radiation is found from the usual scalar radiation transfer equation with the Rayleygh phase 
function and the polarization is calculated as a result of single scattering of known radiation flux before 
the outgoing from the semiinfinite atmosphere. The transformation of this single scattered radiation due 
to Faraday's rotation effect is also taken into account. These formulae give to some extent overestimated 
values for the polarization. The comparison with the exact numerical results for magnetic field directed 
along the normal to the atmosphere shows that for the conservative case (q — dalipa + <^t) = 0) these 
asymptotic formulae give the correct result with the relative error of 10% at (5o = 10. For (5o = 5 the 
error is larger and is equal to k, 20%. For absorbing atmosphere (g ^ 0) these asymptotic formulae are 
less exact then for more important case q — Q. 

The use of these simple asymptotic formulae (see [14]) gives rise to the following expressions for the 
Stokes parameters of the radiation from an optically thick magnetized spherical envelope: 

F, ^ p. (11) 

1 TT 

F ^~9 fj fj /-, 2\ {1 — k^) cos{2if) 
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Here, F is the radiation flux from the 1 cm^ of the envelope surface, a is the outer radius of the envelope, 
R is the distance to an observer, the numerical factors Ji and g depend on q and are calculated in [14]. 
The value 5 cos 9 for the dipole magnetic field is equal to 

5cosQ — (5o [ (3^^^ — 1) cosi?„i + 3^-\/l — pfism-dm cos(/?], (14) 

where Jq — O.SA^M/a'^ is the parameter (2) at the magnetic equator of outer surface of the envelope. 

Further, we shall consider also two cases of the distorted dipole magnetic field: 1) diamagnetic (per- 
fectly conducting) envelope when the normal component of the magnetic field does not penetrate into 
medium, and 2) the envelope with the strong plasma radial outflow when the frozen magnetic fleld 
acquires the radial form (this is the analogy of Parker's spherically symmetric wind for non-rotating 
star). The envelope can be considered as diamagnetic one when the ohmic diffusion characteristic time 
[Tohm ~ ci^ /QDfn with ohmic diffusion coefficient D^ — c^ / 'iircr) is much lesser than others characteristic 
times (the time of turbulent mixing, period of cyclotron rotation and so on (see [15])). 



For diamagnetic envelope the radial component Br of the dipole magnetic field are to be excluded. 
This gives the expression: 

(5 cos 9 = (5o [~ cos -drn (1 ^ M^) + /i-\/l — /i^ sini^m cos ip]. (15) 

For the Parker outflow, on the contrary, only Br must be taken into account: 



(5 cos 9 — 2 5o ip^ cos-Urn + Mvl ~ M^ sa\dm COSip). (16) 

In Figures the cases (14), (15) and (16) are denoted as a), b) and c). 

For small values of the depolarization parameter (5o <C 1 and the most interesting case of conservative 
atmosphere (g = 0) the integrals (12) and (13) can be easily calculated analytically. In this case we have 

Fq (X Sq and Fu (x Sq 't^ Fq. It gives pi ~ \ Fq + Fu/Fi and tan2x — Fu/Fq. Finally we obtain 

PI ^ C 6l sin^ i^m, (17) 

X^ DSf) cosdm, (18) 

where the coefficients C and D for the cases (14), (15) and (16) are equal to: C = 0.765, D = 17.19 for 
the dipole (14); C = 0.0797, D = -51.56 for the diamagnetic envelope (15); C = 0.319, D = 68.75 for 
the radial outfiow (16). It should be noted that X-axis is chosen in the plane "line of sight - magnetic 
dipole" (nM) . The positive x corresponds for an observer to the anticlockwise deviation from this plane 
(nM). If the magnetic dipole M has d = 90°, the sign of the positional angle x is opposite. 

Eqs. (17) and (18) show that pi ex A* and % ex A^. The comparison with numerical calculations 
demonstrates that eqs. (17) and (18) are valid up to Sq = 0.25. 

4 DISCUSSIONS OF THE RESULTS 

Numerical calculations of p; and x from (11)-(13) are presented in Fig. 2-5 as functions of dimensionless 
parameter ^/So = 0.89 A(^m)-\/i3o(G), where Bq = M/a^ is the magnetic field on the surface of optically 
thick spherical envelope. Parameter y/5o changes from up to 10. We calculated polarization for 
q = 0; 0.05; 0.1; 0.5. For values ^/^o > 10 relatively good results can be obtained by extrapolation from 
pi = 100pi{\/So = 10)/So and for positional angles one must be taken x at boundary value v^ = 10 (see 
figures) . 

The peak value of the polarization occurs when M is perpendicular to the line of sight. In this case 
X = 0, i.e. the electric field of the integral polarization oscillates in the plane (nM). For Sq > 100 and 
q ~ we have the asymptotic formula 

Pi = -r 19 

ndo 

where n == 3, 1,2 corresponds to total dipole (14), diamagnetic envelope (15) and Parker's outflow (16), 
respectively. In figures these cases denoted as a), b) and c). 

In the case dm — 90° the expressions (14), (15) and (16) for dcos6 differ only in different coefficients 
before Sq (n ~ 3, 1,2, correspondingly). For this reason, the peak polarization pimax in a-U cases is the 
same but due to the different factors corresponds to different values of ^/5q. For the cases (14), (15) and 
(16), pimax occurs at v^ to be equal to 1.385; 2.4 and 1.895, respectively. Of course, pimax depends on 
the true absorption degree q. The calculations show that pimax increases monotonically from 0.307% at 
g = up to 1.01% at q = 0.5. 

The figures demonstrate that the position of pi max depends weakly on the angle between M and 
line of sight n. For this reason, one can give, the common for all i?„i - angles, approximate formula for 
wavelength Xmax when the peak polarization occurs. For q = the values \/So, mentioned above, give 
rise to the expression: 

2 fif^ 
Xmax ~ -—— inm), (20) 

VnBo 



where n == 3, 1,2 corresponds to the models (14), (15) and (16). Formula (20) allows us to estimate the 
magnetic field Bq on the outer surface of spherical envelope if the value Xmax is found from polarimetric 
observations. For q y^ one can obtain the analogous formula. 

The procedure is the following: from b) - curve one can find the value (\/5o)max, divide it by VO.S == 
0.894 and take the number instead 2.86 in (20). For q = 0.5 such procedure gives 2.07 instead of 
2.68, i.e. the existence of absorption changes the estimation of Bq only slightly. Roughly, one can take 

1/2 

Xmaxi.tJ'm) « Bq (G). 

Let us discuss now how to know from the polarimetric data whether the envelope is optically thin 
or optically thick. First of all, remember that the polarization degree of radiation, scattered in optically 
thin envelope, as a rule is sufficiently higher than that for optically thick envelope. But, the dilution 
effect from the nonpolarized radiation can decrease this difference. For example, the question arises from 
the analysis of AGN polarimetric data how the dilution effect from the neighboring stars is significant. 
It seems more informative to consider the spectral differences of polarization in these cases. The optical 
depth of the envelope, in principle, depends on wavelength. It can give the jump of positional angle 
at some A. This means the transition from optically thin envelope to the optically thick one (see the 
monograph [1]). The calculations in [8,9,11] for optically thin envelopes are made only for pure dipole 
field (the case (a) in our figures). Therefore, we shall compare only upper figures with the curves in Fig.l. 

The curves of the polarization degree pi in both types of envelopes are qualitatively similar. First, 
the increasing of pi ex A^ takes place, then the peak exists and finally polarization degree decreases. 
This decreasing depends on the explicit form of number density Ne{r). For envelopes with permanent 
concentration [Ng = const) we have pi ex A^ at large values of A. This dependence coincides with 
the corresponding asymptotic behavior in the case of optically thick envelope. But, most probably the 
number density of free electrons falls with the removal from the central source. If N^ (r) is proportional 
to r~^, then in the case of optically thin envelope the polarization degree decreases as ex l/vA whereas 
optically thick envelope gives pi ex 1/A^. 

Thus, the form of the polarization spectra at large wavelengths allows to say whether the envelope 
is optically thin or thick. Unfortunately, the value of polarization at large wavelength can be very small 
(see (19)) and be undetectable from observations. Besides, the asymptotic region of large wavelengths 
depends on unknown magnetic field Bq on the outer surface of the envelope. 

To estimate the magnetic field Bq for both types of envelopes we need to know the spectrum of 
polarization near the pimax value, i.e. we need to know Xmax- After that we can use the relation (20) in 
the case of optically thick envelope. Analogous expressions for optically thin envelopes are presented in 
[8,9]. For most probable distribution N^. ex r~^ the spectrum of polarization of optically thin envelope 
has the peak at the following Xmax'- 

^ 4.65 ^ 4.65 rV^ 

'^max ~ r-j= r-j= • 3/2 • l^-Lj 

V i>s Tenv V -DO Tenv ^ ' 

Here, Bg — M/R^ is the dipole magnetic field at the surface of central radiation source. The dipole law 
gives Bq = BgR^/a^- Because of i?s < a and Tenv < 1, the estimation (21) gives rise to larger value 
for Xmax, than the estimation (20). It should be noted that for very large Bq the value Xmax can be in 
the X-ray range of the wavelengths. If the value Bq is assumed to be known then the estimations from 
asymptotic formulae can be made without knowledge of Xmax - 

The spectra of the positional angles x for optically thin and thick envelopes are different. Generally 
speaking, this fact gives an additional opportunity to distinct these both cases. However, if the observed 
polarization corresponds to large wavelengths (when Sq ^ 1) then in both cases the angle x has practically 
permanent value and sign. In this case we can say nothing to resolve the dilemma. Really, the spectrum 
of X is useful for this aim, when it is known near Xmax, where the behavior of spectra for both cases is 
very different. 

In the region oi Sq <^ 1 the x-spectrum is practically useless as it seems that the signs of x ^^c 
different in these cases. The sign of x depends on the dipole orientation which, as a rule, is unknown. If 
M changes orientation to opposite one, then x changes its sign. However, if from some reasons we know 
the M - orientation, then, in contrary, the observed sign of x demonstrate directly what type of envelope 
exists - optically thin or optically thick. 



Table 3: The net polarization pi of hot stars for various magnitudes of equatorial magnetic field. 
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5 SOME ASTROPHYSICAL APPLICATIONS 



5.1 



Polarization of radiation from OB and WR stars: 
magnetic field 



the estimation of a 



The hot stars are characterized by very strong, ~ M > 10^ M©, outflow of a matter in form of the dense 
stellar wind (sec [16-18]). The Thomson optical depth t of such wind can be both less {( Puppis, Dencb) 
and greater (of the order) than unity (P Cyg, WR 40) (see [17]). The problem of the determination of 
surface magnetic field for these objects exists for a long time. In the paper [19] the positive answer is 
given for the question whether dipole magnetic field can exist at the surface of these stars. 

The attempts to measure the magnetic fields of bright stars by the traditional circular polarization 
technique were unsuccessful [16,20]. Particulary, in the paper [20] the upper limit of the magnetic field 
component along the line of sight was obtained for the star 04I(n)f ( Puppis: _B|| < 2006*. The sufficient 
variability of the Stokes parameter V was not also found. One can expect that future Spcctro-Polarimetric 
Interferometer (SPIN) will be useful. 

Using our calculations (Fig.l and 2) one can estimate the polarization in V-band for a number of stars 
from the list of objects of the project SPIN [17]. So, assuming that equatorial magnetic field Be ~ lO^G 
we obtain Sq w O.SXyBe « 24.2. If Be « 10 then we have Sq « 2.42. The values of expected polarization 
degree are presented in Table 3 (for the inclination angle i — 90°). In estimations we have chosen the 
wind density distribution from the paper [21]. 



5.2 Polarization of radiation from the system Cyg X-l/HDE 226868 

In 1975 Nolt et al. [22] have discovered the intrinsic linear polarization of the optical radiation of 
close binary system Cyg X-l/HDE 226868 with a compact component being a black hole with the mass 
~ IOMq. The polarization has variations with the amplitude ~ 0.2%. The time behavior of polarization 
is rather complex. Besides the period of the binary system P = 5''.6, the polarization varies with 39^^ and 
78'* characteristic times [22]. The possible mechanisms of the polarization of this system were discussed 
in [26]. 

If the polarized radiation arises in the accreting matter around the black hole the value of intrinsic 
polarization is to be some tens times higher than the observed value. It seems there exists very strong 
dilution by the unpolarized emission of supergiant HDE 226868. The luminosity of this star is estimated 



as Ln 



(1^3) X 10^9 erg fi-i 



whereas the X-ray luminosity of Cyg X-1 is about Lx < 8 x 10 ergs 



[23] . If optical radiation of the accretion disk arises as a result of the transformation of X-ray emission 
of the black hole then the value of intrinsic polarization have to be near the value pi < 10%. Optically 
thick accretion disk with most probable inclination i = 30° can not produce such high polarization. For 
this reason if the polarization originates into accreting matter, then this polarization can arise only as a 
result of the scattering in optically thin coronae or in some sort of outflow from the accretion disk (wind) . 
At the same time, it is not possible to exclude the mechanism of the production of small (at the level of 
some tens percents) intrinsic polarization from the star HDE 226868 in a magnetized stellar wind. 

We consider both mechanisms and the corresponding consequences. 

In the case of optically thin magnetized corona around the accretion disk the large, ^ pi ^ 10%, 
polarization can be produced if the central source is point-like or the coronae itself is very vast. For such 
model one can use the calculations of Dolginov et al. [1]. According to their calculations pi ^ 10% can 
be obtained at SoTenv ~ 10 and when magnetic field is perpendicular to the line of sight. For i = 30° 



it means that the magnetic field Ues almost in the plane of the disk. If one supposes re„„ « 0.1 then 
the condition Sq w 100 (remember that 6oTenv ~ 10) gives rise to the estimation B^ w 500 G. In a 
number of models the inner radius of a hot coronae is determined as ^ 100 i?g, where Rg ~ 25M/c^ is 
the gravitational radius. In this case for dipole magnetic field follows the estimation B{3Rg) ^ 10^ G. 

Let us now consider the situation when the linear polarization is generated in the extended stellar 
wind of the optical component HDE 226868. The existence of such wind was confirmed by observation 
of modulated radio emission from Cyg X-1. The radio emission of the relativistic jet of the black hole 
has variable absorption in the plasma of the stellar wind. Modulation arises in the process of the orbital 
motion of the black hole [27] . 

Supposing that in stellar wind N^ ^ r^^ and M « 2 x lO^^MQyear^^, Eg = IORq and T4o ~ 
1850A;tos^^ [27], one can estimate the Thomson optical depth of the wind as tt ~ 0.1. According to 
Fig.l the polarization pi « 0.2% corresponds to the angle between magnetic dipole M and line of sight 
being ^ 65° and the parameter \^5st^ ~ 3. It gives the estimation Bs ~ 350 G. Such value of extremal 
magnetic field for the OV supergiant seems to be rather large. It is interesting that in this case the dipole 
axis lies in the orbital plane and one has not any problems to explain the variability of the polarization. 

5.3 Intrinsic polarization of SS 433 

The famous galactic X-ray binary system SS 433 is the microquasar with two relativistic {V ~ 0.26c) 
strongly collimated (i? w 1°) jets. The system is characterized by 3 periodic motions: the orbital motion 
with the period Po = 13^^.082, the precession with Pp = 162''. 5 and nutation with P„ = 6^^.28 [23,28]. 
Recent spectral observations [20] confirm the existence of spectral absorption lines of the A7Ib supergiant. 
The mass ratio of the relativistic component and an optical star is about q = Mx/My — 0.57±0.11. The 
masses of every component are estimated as My = (19 ± 7)Mq, Mx = (11 ± 5)Mq [28]. The kinematics 
of relativistic jets allows to estimate the inclination angle of the orbital plane i = 78°. 82 ± 0°.ll. The 
precession angle of a disk is equal to 9 = 19°. 80 ± 0°.18. 

The variable, i.e. intrinsic, polarization of SS 433 first was observed by McLean and Tapia [30], and 
later, in B, V, R, I bands, by the observers in KrAO [31]. In 1993 the polarimetric observations of 
SS433 were obtained in UV range by high speed photometer-polarimeter (HSP) from the cosmic Hubble 
telescope [32]. The polarization degree in UV range is rather high pi = (13.4 ± 4)% and it demonstrates 
the time variability. Such a high degree of polarization can not originate in an accretion disk. Most 
probably it is a result of the scattering of radiation by free electrons in a plasma outflow (for example, in 
a hot corona or wind). An example of such plasma can be, also, optically thin advective accretion flux 
around the black hole [33]. 

We calculated the expected polarization degree of the radiation scattered in a spherical stellar wind 
with the Parker's type magnetic held 

Br = BJRA/rf cos 6*; B^ = Bt— sin 61; Bg^Q (22) 

r 

where B^ is the magnetic field at magnetosphere's pole, Bt is the value of toroidal component, Ra is the 
radius of the magnetosphere. There was used the technique described in the papers [1,8-10]. 

The results are presented in Fig. 6. They are compared with the observational data of Dolan et al. 
[32]. We assumed that the Thomson optical depth of the coronae is near unity, i.e. tt < 1- 

The results of UV polarimetric observations corresponds rather well to the theoretical curve if one 
assumed Bp = IG and Bt = 600 G. The disagreement between observational polarization and the theory 
in optical range can be explained by the existence of the dilution from intrinsic non-polarized radiation 
of the accretion disk. However, the Thomson polarization of radiation outgoing from plane-parallel disk 
at i w 80° is rather high {pi « 5%, [13]), the strong {Bt « 600 G) toroidal magnetic field depolarizes 
fully the scattered radiation because of the depolarization parameter do is very large {Sq > 500(A/1/zto)^, 
see (21)). 

It is interesting that our estimation of the magnetic field near the magnetosphere allows us to obtain 
the additional information about the parameters of the accretion to the black hole itself. 

The modern theory of rotating black holes determines the magnetosphere as a region from which the 
relativistic jets outflow occurs and, simultaneously, arises the hot coronae [4,34]. According to [34], the 
radius Ra ~ 5 x 10^i?g in the case of pure thermal coronae and Ra ~ 300i?g for the, so called, hybrid 
coronae. As a result one can estimate the value of magnetic field in the vicinity of black hole, near 
the last stable orbit 3Rg {Rg is the gravitational radius). If one assumed that the estimation (21) for 



Parker's type of magnetic field is valid up to 3Rg, then we obtain the estimation of the magnetic field of 
the black hole in SS 433 as Bmax ~ 10'* -^ 10® G, depending on the type of coronae. This value is in a 
good agreement with that in the model of Blandford and Znajek [4]. 

It should be noted that if one takes Ra ~ 10*° cm (as in [34]) then the value of magnetic moment fj, 
of the compact component in SS 433 can be estimated as /i w 10'^^ Gcm^. This value is near the value 
of /i for another microquasar GRS1915+05 obtained by Robertson and Leiter [2] from very different 
consideration. 

Finally, if, according to [33], one considers that the magnetic force lines strongly follow the motion 
of the plasma in the accretion disk, then we can estimate the known viscosity factor a of the Sunyaev- 
Shakura theory [35] . According to 

— = — K a ('231 

where Ur and U^ are the radial and Keplerian velocities in accretion disk, correspondingly. Our estimation 

■"v 
[36]. 



of Br and B^p gives rise to a « 2 x 10 '^ what is practically coincides with the theoretical estimation in 



5.4 Polarization of radiation of supernovae 

Recently a lot of observational data has been obtained giving an evidence that the radiation from super- 
novae has a sufhcient polarization (see the review [37]). In particular, the spectropolarimetric observations 
of young supernovae, obtained at Kek-telescope, demonstrate that the intrinsic polarization exists for all 
the types of supernovae [38]. These observations, for example, show that the supernovae 1997dt (la), 
1998T and 1997dq (both of the type lb), 1997ef (Ib/c - pecular), 1997eq (Iln), 1997ds (II-p) have intrinsic 
polarization. It was shown [39] that the supernovae of types II and Ib/c are polarized at the level p; ss 1%. 
Some of them have higher polarization. It is interesting (see [39]) that the lesser is the hydrogen envelope 
and the deeper layers are seen the larger is the value of polarization. Such dependence was analyzed by 
Hofflich et al. [40]. It proved to be the case that the observed polarization is higher for SN II as compared 
with SN Ic. 

It is natural, that many authors explain the observed polarization as due to the scattering of radiation 
in asymmetric envelope of the expulsed matter (see, for example, [37]). 

There exists another physical mechanism that can be responsible for the noticeable net polarization 
of SN radiation without asymmetric distribution of the exploded matter. 

Such mechanism is the considered above integral effect of the Faraday rotation of the polarization plane 
of the radiation scattered in the spherically symmetric envelope of exploded matter (blow shock) with 
a magnetic field. Namely this effect produces the intrinsic polarization. The measured net polarization 
allows to estimate the magnetic field strength in the region of propagating blow shock. 

Now let consider some examples. The supernova SN Up 1999em was observed at epochs of 7, 40, 49, 
159 and 163 days after the explosion and its radiation was to be polarized at these epochs (Leonard et al., 
2001). The detected net polarization was pi ^ (0.2 4- 0.5)%. This magnitude is quite well corresponding 
to the results of our calculations presented at Fig. 2, the polarization maximum corresponding to the 
Faraday depolarization parameter Sq ^ 4. For the visual magnitude (V band - Xeff = 0.55fnn) this value 
of the depolarization parameter allows to determine the magnetic field strength as _B ss 16.5 G. 

Another example is SN la. The typical value of the intrinsic polarization of such kind supernovae is 
to be pi « 0.7% (Kasen et al., 2003). This value of polarization is better corresponding to our results 
presented at Fig.l, i.e. to the case of scattering into the magnetized optically thin envelope. In this case 
the real depolarization parameter is (Sstt)^''^ , where ry <C 1 is the optical thickness of an envelope with 
respect to the electron scattering. The polarization magnitude p; ~ 1% corresponds to the depolarization 
parameter value of (dsTr) ^ 1- In this case the magnetic field strength is Bs ^ 4/tt G. For tt ~ 0.1 the 
magnetic field strength in the shock region can reach the strength ^ 40 G and even ^ 100 G. 

The calculation of non-sphericity of an envelope is not principal problem. Some cases of generation 
of polarized radiation in a magnetized non-spherical envelope was considered in [1,14]. 
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6 CONCLUSIONS AND OUTLOOK 

We represented the results of calculations of the integral polarization of radiation emerging from the 
optically thick spherically symmetric envelope with the dipolc magnetic field of a central radiation source 
(a neutron star; extended supernova envelopes produced by the supernova explosion at the initial stage 
of expansion; accretion flows around quasars and active galaxy nuclei). We considered also the cases 
of diamagnetic (perfectly conducting) plasma envelope when the normal component of magnetic field 
does not penetrate inside a disk and radially evaporated envelope. The comparison of the results of 
polarization between the various types of envelopes have been made. 

The results of our theoretical calculations can be used for the analysis of polarimetric observations of 
hot stars, including WR stars, X-ray binaries (Cyg X-l/HDE 226868, SS 433), supernovae and etc. We 
plan to use the results of our calculations to an analysis of polarimetric data of quasars and AGNs (see, 
for example, [44]). 

This work is partially supported by the RFBR grant 03-02-17223, the Program of the Prezidium of 
RAN " Nonstationary Phenomenae in Astronomy" , the Program of the Department of Physical Sciences of 
RAN "The Extended Structures..." and by the Program of Astronomy of Russian Science and Education 
Ministry. 
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Figure 1: The polarization degree and the positional angle of the integral radiation of optically thin 
spherical envelope in the magnetic dipole field. The upper panels corresponds to the case Ne ^ r^^ and 
the low ones to the envelope with constant electron number density Ne = const with 77 = 5. The numbers 
near the curves denote the angle (in degrees) between the magnetic dipole M and line of sight n. 
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Figure 2: The polarization degree and the positional angle of the integral radiation of the optically thick 
magnetized for the conservative atmosphere (q — 0). The numbers denote the angle 'dm (in degrees) 
between the magnetic dipole M and the line of sight n. The cases a) , b) and c) correspond to the total 
dipole field, diamagnetic atmosphere {Br = 0) and Parker's outflow field (Br j^ 0), correspondingly. 
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Figure 3: The same as in Fig. 2, but for q — 0.05. 
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Figure 4: The same as in Fig. 2, but for q = 0.1. 
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Figure 5: The same as in Fig. 2, but for q = 0.5. 
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Figure 6: The spectrum of integral linear polarization of X-ray binary system SS433: comparison of the 
observational with the results of theoretical model (solid curve) at tt ~ 1.0, Bpoi ~ IG, Bror = 600G. 
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